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ABSTRACT 

We present the photometric and spectroscopic study of the very luminous 
Type Iln SN 2006gy for a time period spanning more than one year. The evo- 
lution of multiband hght curves, the pseudo-bolometric (BVRI) hght curve and 
an extended spectral sequence are used to derive constraints on the origin and 
evolution of the SN. 

A broad, bright (Mr ~ —21.7) peak characterizes all monochromatic light 
curves. Afterwards, a rapid luminosity fading {•jr ~ 3.2 mag (100 d)"^) is fol- 
lowed by a phase of slow luminosity decline (7^^ ~ 0.4 mag (100 d)^^) between 
day ~170 and ~237. At late phases (> 237 days), because of the large luminosity 
drop (> 3 mag), only upper visibility limits are obtained in the B, R and I bands. 
In the near-infrared, two K-band detections on days 411 and 510 open new issues 
about dust formation or IR echoes scenarios. 

At all epochs the spectra are characterized by the absence of broad P-Cygni 
profiles and a multicomponent Ha profile, which are the typical signatures of 
type Iln SNe. Ha velocities of FWHM ^ 3200 km s"^ and FHWM ^ 9000 km s"^ 
are measured around maximum phase for the intermediate and high velocity 
components, respectively, and they evolve slowly with time. 

After maximum, spectroscopic and photometric similarities are found between 
SN 2006gy and bright, interaction-dominated SNe (e.g. SN 1997cy, SN 1999E 
and SN2002ic). This suggests that ejecta-CSM interaction plays a key role in 
SN2006gy about 6 to 8 months after maximum, sustaining the late-time-light 
curve. Alternatively, the late luminosity may be related to the radioactive decay 
of ~ 3M0 of s^Ni. 

Models of the light curve in the first 170 days suggest that the progenitor was 
a compact star (R ~ 6 — 8 • lO^^cm, Mgj ~ 5 — I4M0), and that the SN ejecta 
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collided with massive (6 — IOM0), opaque clumps of previously ejected material. 
These clumps do not completely obscure the SN photosphere, so that at its peak 
the luminosity is due both to the decay of ^^Ni and to interaction with CSM. 
A supermassive star is not required to explain the observational data, nor is an 
extra-ordinarily large explosion energy. 

Subject headings: supernovae: individual (SN2006gy), stars: circumstellar mat- 
ter, mass loss — techniques: photometric, techniques: spectroscopic 



Introduction 



Textbook stellar evolution theory explaii is that Type II SNe are produced by t he core 
collapse of H-rich stars with masses >8M0 JSranch et al.lll99ol : krnett et al.lll996h . Con- 
firmation of this scenario comes from the possible identification of the progenitors of a 
few SNe II, includ i ng SN 1987A (lArnett et al.l Il989l ) and several more recent events (e.g. 
Smartt et all 120041 : 1 Van Dvk et al.l[2003h . 



On the other hand, a fully consistent picture of a SN progenitor evolution and explosion 
is still missing. Parameters like progenitor radius, ejecta mass, explosion energy, asymmetries 
and radioactive elements yield, all contribute to determining the SN display. 

One of the most uncertain ingredients is the progenitor mass-loss history. Although 
some constraints on the progenitor mass loss can be derived from models of the SN light 
curves and spectra, direct measurements are only possible when the circumstellar material 
(CSM) becomes visible. 



Denser CSMs and hence higher mass loss rates (~ 10 ^ Mq yr ^. ISalamanca et al.ll2002l ) 
are required to explain the su dden halt in the la te-time lu minosity decline of s everal type 
II Linear SNe (e.g., SN 1979C, lBraiidi_elalJll98ll; SN 1980K. iMontes et al.lEggsl : SN1986E, 
Cappellaro et al.lll995l ). 



Mass loss plays a key role in type Iln SNe (ISchlegel et al.lll990l ). The spectra of this 
class of SNe are characterized by emission lines with multiple components, which may range 
from very high (~ 20000 km s^^) to low velocities (a few hundred kms ~^) with no associated 
(broad) P-Cygni absorptions. Such events may be very energetic (e.g.. lAretxaga et al.lll999l ) 
and are characterized by a slow luminosity evolution beginning soon after discovery. It 
is generally believed that the shock produced when the high velocity ejecta impact on a 
relatively dense CSM causes the conversion of part of the ejecta kinetic energy into radiation. 
Depending on the CSM density distribution, the strong CSM-ejecta interaction may last for 
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months (e. fi., SN1994W, 
SN1995N, 



Fransson et al 



Chugai et alJ 120041) to years fe .g.. SN198 8Z. IXuratto et a 



2002 and lZampieri et al.ll2005l : SN 1995G. IPastorello et al. 



1993 



2002|). 



The recent peculiar SN Il n 2006gy event attracted much interest. Discovered on 2006, 
September 18th in NGC 1260 JOuimbvl 120061 ) . SN2006gy was initiallv classified as a SNII 
( iHarutyunyan et al.ll2006l ) and shortly thereafter SNIIn JPolev et allbood l Although 
its spectroscopic features were not unprecedented, the photometric behavior was. A very 
bright luminosity peak (Mr ~ —22, ISmith et al.l 120071 . hereafter S07) was in fact reached 
only ~ 70 days after explosion (estimated by backward extrapolation of the light curve, cr. 
§2). The amount of energy radiated in visible light during the first 200 days (> lO^^ergs"^, 
S07) was larger than in any previously observed SN, either core-collapse or thermonuclear. 



On the other hand, in contrast with other bright SNe Iln (jZampieri et al.ll2005l . I Chandra et al 



20051 ). only a weak and soft X-ray emission was detected by Swift and Chandra near the 
epoch of optical maximum (S07). The absorption-corrected, absolu te X-ray luminos ity in 
the band 0.65 - 2keV was 1.65 x lO^^ergs-^ (S07). It was argued Jofek et al.l[2007l . S07) 
that if direct radiation from the ejecta-CSM interaction was the cause of the extraordinary 
SN luminosity, the X-ray fiux should have been several orders of magnitude larger. A pos- 
sible explanation is that the interaction region might have been hidden because of a very 
large optical depth. Alternatively, the source for the optical luminosity should be searched 
elsewhere. 

Actually, the first model proposed by S07 was the explosion of a very massive star 
(> IOOMq) via pair-instability phenomena. Such a violent explosion would cause the com- 
plete disintegration of the core and the ejection of a huge amount of ^^Ni, possibly up to 
22M0. In such a scenario the energy input from the radioactive decay chain ^^Ni -^^^ Co — Fe 
accounts for the observed luminosity. 



However, late luminosity measurements presented by lSmith et al.l (120081 ) (hereafter S08) 
and in this work set a much estimate to the possible ^^Ni mass. Other models call for 
the conversion of some of the ejecta kinetic energy into radiation , via a collision with a 
massive (~ IOMq) and highly opaque (r ~300) circumstellar shell (jSmith fc Mc Crayl 12007 
and S08). This scenario provides no direct information on the real nature o f the explosion 
which, in principle, could even be a thermonuclear runaway (jOfek et al.ll2007l ). However, this 
is difficult to reconcile with the presence of the massive circumstellar shell. Thus, the core 
collapse of a massive star is still the most appealing scenario. Alternatively, we may have 
witnessed the collision between high velocity shells originated in subsequent outbursts of a 
very massive st ar undergoing struct ural instabilities caused by pair production {pulsational 
pair-instability, IWoosley et al.l 120071 ). In any case, the presence of massive shells (~ IOM0) 
at a large radius (i.e., a few 10^^ cm) is required to explain the observed high luminosity. 



- 5 - 



Such a large shell mass, couple d with a high density, may result in long photon diffusion 



times (jSmith fc Mc Crayl 120071 ). which could explain the broad peak of the light curve of 



SN2006gy. 

These scenarios can explain the light curve evolution during the first 5 months. However, 
the presence of a diffusion process implies fairly rapid fading of the luminosity after maximum 
(i.e., a time of the order of the photon diffusion time-scale). Observations instead indicate 
that the luminosity does not decline rapidly after maximum. On the contrary, the light 
curve decline slows down between day 170 and day 237 (see §3). This requires an additional 
energy source. Again, r adioactive material has been proposed (~ 8Mq of ^^Ni according to 



Smith fc Mc Crayl 120071 ). but new doubts have arisen after SOS reported a late-time drop of 
the optical luminosity and two infrared SN detections, which may support the formation of 
dust in the ejecta or the presence of IR echoes. 

In this work we present new data, which include observations at all relevant phases, 
from discovery to more than 1 year later. Multiband light curves and an extended spectral 
sequence are shown and discussed. By comparison with other SNe and by means of modeling 
we try to verify to which extent this SN is really extra-ordinary and to provide new constraints 
for the progenitor and the explosion. 



2. Observations 

Optical (BVRI) and near-infrared (JHK') images of SN2006gy were acquired at TNG, 
NOT (La Palma, Spain) and the Copernico 1.82m Telescope on Mt. Ekar (Asiago, Italy) 
over a period spanning more than 500 days from discovery. Optical spectroscopy was also 
performed, up to ~389 days (see Table 1 for a complete log of the observations). 

Since the SN is located very close to the nucleus of the host galaxy, template subtraction 
was required for photometry. We used archival B, R and I band images of NGC1260 acquired 
at the Jakobus Kapteyn Telescop^ (La Palma, Spain), and a V band image taken at the 
Schmidt Telescope at Kiso Observatorj]^ (Japan). Additional information about the template 
images is given in Table 2. 

All images were de-biased and fiat-field-corrected. A local sequence of stars in the 
SN field was calibrated using observations of standard stars obtained during photometric 



^web archive http:/ /archive. ast.cam.ac.uk/ingarch/ingarchold. html 
^web archive http://www.ioa.s. u-tokyo.ac.jp/kisohp/ 
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nights. Template subtraction was performed using ISIS (lAlard &: Luptonlll998l ). and the SN 
magnitudes were measured on the subtracted image with a Point Spread Function-fitting 
technique, using custom-made DAOPHOT-based routines that were adapted specifically for 
supernova photometry. 

For spectroscopy, all scientific exposures were acquired at low airmass and positioning 
the slit along the parallactic angle. Wavelength calibration was accomplished with arc-lamp 
exposures and checked against the night-sky lines. The flux was calibrated using instrumental 
sensitivity functions obtained from observations of spectrophotometric standard stars. These 
were used also to remove telluric absorptions from the spectra. In order to improve the signal- 
to-noise ratio, separate spectra taken during the same nights were combined. Finally, flux 
calibration was checked against photometry. If necessary, a constant multiplicative factor 
was applied to correct for flux losses caused by slit miscentering or non-photometric sky 
conditions. 

The spectrum acquired on 2006 December 19th with the Ekar 1.82m telescope required 
further adjustments because of the poor seeing conditions and the residual contamination 
from the galaxy background. The latter was removed using the spectrum at 389 days, where 
the SN is not detected, ClS db background template. 

Throughout this paper, for the sake of simplicity, phase refers to the same reference 
epoch as in S07, JD=2453967 (2006 August 19.5 UT). This was derived from a backward 
extrapolation of the rising branch of the light curve. S07 refer to this date as the explosion 
epoch, but this term may be misleading in view of some of the proposed scenario^. 



3. Photometry 



To compute the absolute magnitudes of SN 2006gy some assumptions regarding the host 
galaxy distance and SN extinction are needed. 

Lacking other indicators, the distance to NGC 1260 was estimated from the Hubble law, 
Vrec/Ho, where Vrec = 5822kms~^ is the host galaxy recession velo city corrected for Virg o 



cluster infall (from HyperLEDM) and Hq = 72 ± 8 kms^^Mpc"^ dPreedman eraPboOlh . 
These values imply a distance modulus /x = 34.53, equivalent to a distance of 80.86 Mpc. 



•^In the shocked-shell diffusion model ijSmith &: Mc Cravll2007l ) the detection of the SN emission occurs 
a few weeks after the real explosion. In the scenario suggested bv lWooslev et al.l (|2007r ) there is not even a 
SN explosion, but only an outburst release of matter. 



*http:/ /leda. univ-lyonl.fr/ 
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We adopted a G alactic extinction towards NGC 1260 of E(B — V)gai = 0.16 (AB,gai = 0.69, 
Schlegel et al.lll998l ). An estimate of the extinction in the host galaxy was obtained com- 
paring the spectra of SN2006gy to those of SNII 2007bw, another peculiar and bright 
SNIIn, photometrically similar to SN2006gy. This yields E(B — V)host — 0.4, assuming for 
SN2007bw little or no extinction. It is interesting to note that on the 42 days spectrum the 
EWs of NalD due to the galactic and interstellar absorption are 2. 2 A and 5.5A, respectively. 
Assuming for NGC 1260 a gas-to-dust r a tio al ong the line of sight as in our Galaxy and 
adopting the extinction of ISchlegel et al.l (Il998l ). the derived internal absorption is fully in 
agreement with that derived by comparison with SN 2007bw. 

Therefore the total color excess is Etot(B — V) ~ 0.56. This is comparable to the esti- 
mate of S07, i.e. Etot(B — V) ~ 0.48, which wa s also obtained by comparison with SNelln, 
and slightly smaller than the value adopted by lOfek et al.l ( 120071 ) . Etot(B — V) ^ 0.7. 



The BVRI absolute light curves of SN2006gy are plotted in Figure 1. In all bands 
the light curve exhibits a slow increase to maximum, which is reached ~ 70 days after the 
reference epoch. The peak magnitudes are B ~ —21, V ~ —21.4, R ~ —21.7 and I ~ —21.5. 
Such an extended, plateau-li ke peak was noted for a type Iln SN only in the case of SN 
2005kd JTsvetkov et al.lboosh . 



Between day ~ 100 and day ~ 170, the light curve declined relatively rapidly (75 ~ 
S.Omag (100 d)~^, 71/ ~ 3.1, 7i{ ~ 3.2, 7/ ~ 2.8). Then, from day ~ 170 onwards the 
light curve evolution suddenly flattened: in the following ~ 70 days the decline was only 
7ij ~ 0.4mag (100 d)~^. When the SN could be observed again after solar occultation, its 
luminosity was below the detection limit in the optical bands. A limit was obtained placing 
artificial stars of different magnitudes at the SN position. Despite the long exposure times, 
only relatively bright upper limits were derived because the SN is very close to the nucleus of 
the galaxy (cf. S07, Figure 1). The derived apparent magnitude hmit on day 389 is > 20.3 in 
R. Optical upper limits were obtained also at 423 days, when we derive B > 21.0, R > 21.5 
and I > 19.75. These measurements imply a new steepening in the luminosity decline after 
day 237. 



Guided by the evolution of other SNe (e.g., SN1998S, 


Pozzo et al. 


200 


4|or S 


Tominaea et al. 


2007; 


Smith et al. 


2008; 


Mattila et al. 


2008 


Di Carlo et al. 


bood) 



we con- 



sidered the possibility that at late epochs a significant fraction of the bolometric luminosity 
could be emitted at IR wavelengths. To test whether this was the case, late observations of 
SN2006gy were obtained with NIGS at the TNG, on 2007, October 5 (JHK' bands) and on 
2008, January 12 (K' band only). We could not apply the template subtraction technique 
in the near-infrared because the available pre-discovery images of the host galaxy retrieved 
from the 2MASS archive are not deep enough to be compared to the TNG images. There- 
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fore, we had to rely on the PSF fitting technique which, given the SN position, has a large 
uncertainty. Photometric cahbration was performed adopting field star magnitudes as listed 
in the 2MASS Point Source Catalogue!. 

The SN was not detected in the J and H bands, for which we could only estimate upper 
limits, J > 17.0 and H > 16.5. Instead, a point source was detected in the K' band at the SN 
position (Figure 2). The SN was measured at K ~ 16.0 ±0.5 on day 411 and K ~ 16.3 ±0.5 
on day 510. These values are ~ 1 mag fainter than those measured by SOS at similar epochs 
(K=15.1 ± 0.1 and K=15.4 ± 0.1 on day 405 and 468, respectively). Even allowing for the 
large error bars, the two sets of measurements do not agree, probably because of a different 
calibration. 



3.1. Infrared emission and bolometric luminosity 



Given the K'-band detection of the SN at 411 and 510 days, it cannot be excluded that 
at late phases a considerable amount of flux is emitted in the near-infrared. 

S08 suggested two possible sources of the late K'-band luminosity. The K'-band emission 
could be associated with an IR light echo from circumstellar dust, for which the input energy 
is the light emitted by the SN near maximum. In this case the IR flux should not be 
considered when computing the late-time bolometric light curve. Alternatively, the IR flux 
may originate from circumstellar dust heated by an instantane ous energy supply ( radioactive 
decay or on-going CSM-ejecta interaction), as was suggested by lPozzo et al.l (120041 ) to explain 
the late phase photometric data of SN1998S. 

In order to get some constraints on the total emission from dust at ~ 411-423 days we 
assumed a black bo dy energy distr ibution multiplied to a factor 1/A, as an approximation 
of what reported in ISpitzer I (119981 ). and normalized it to the observed K-band flux. Given 
that we have no constraints on the dust temperature, we adopted th ree values includiri g 
T=1200K, the dust temperature in the ejecta of SN 1998S derived by IPozzo et al.l (120041 ). 
For each value we plotted the spectral energy distribution (SED) of the associated emission 
(Figure 3) and integrated over the entire wavelength range from Ak = 2.16 fim to A = oo. 

It is interesting to note that the K-band luminosity of SN 1998S measured at similar 
epochs (K=13.8 at ~464 days) would differ from that of SN 2006gy by a factor 1-5 (adopting 
for SN2006gy K=16.3 from this work or K=15.4 from S08) if scaled at the same distance. 
Therefore, given that the two fluxes are of the same order, it is plausible that any mechanism 



''http: / /tdc-www. harvard.edu/software/catalogs/tmpsc. html 
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explaining the IR emission of SN 1998S can work also for SN2006gy. 

From our multiwavelength photometry we can derive the pseudo-bolometric luminos- 
ity evolution of SN2006gy integrating the flux in the optical bands (BVRI). The pseudo- 
bolometric light curve is shown in Figure 4, compared to those of the type II SNe 1987A, 
1995G, 1997cy, 1999E and 2005gj. The pseudo-bolometric luminosities which include the 
dust contribution in the near IR are represented with plus symbols. It is remarkable that, 
at about 6 to 8 months, the luminosity and decay rate of SN 2006gy become comparable to 
those of other events, in particular to SNe 1997cy and SN 1999E. We will come back to this 
issue in 54 and 55.2. 



4. Spectroscopy 



The spectral evolution of SN 2006gy from day 37 to day 389 is presented in Figure 5. 
The spectra have been de-redshifted and corrected for extinction. 

The early spectra (days 37 and 42) show the typical features of SNe Iln, namely a blue 
continuum and strong Ha and H/5 emission lines, without broad P-Cygni absorptions. 



Using GELATO (iHarutyunyan et al.ll2008l ). the automatical spectra comparison tool 
applied to the Asiago Supernova Archive (ASA), the best match for the earliest avail- 
able SN2006gy spect rum (phase 37) is found with SN 1995G (~36 days since discovery, 
Pastorello et al.ll2002l ). which is generally considered a typical SNIIn, although there are 
differences in the strength of the lines emission. 

Three months after discovery, the spectrum of SN2006gy became similar to those of 
other well-studied SNe II. In the high resolution spectrum taken on day 96 by S07 (their 
Figure 4) narrow absorption lines of Fe II (multiplets 42, 48, and 49 at 5000-5400 A and 
multiplet 74 in the region 6100 -6500 A) are evident. A similar narrow line forest was identified 



in the spec tra of SN 1999el flPi Carlo et al.lb ooj), SN 1995G JPastorello et all l2002h and 



SN 1994W ( SoUerman et al. 19981 : Chugai et al. 2004). In all cases these lines are associated 



with slowly expanding, unperturbed material surrounding the star. 

Despite the lower resolution and S/N ratio, our spectrum at phase 122 days is broadly 
consistent with the features discussed by S08. At this phase the Ha flux has decreased by a 
factor of 2 with respect to the first spectrum, while H/? emission almost disappeared. 



At a phase of 174 days the near-IR Ca II triplet is strong in emission. For this epoch a 
good spectral match is obtained with SNel997cy, 1999E (Figure 6), 2002ic and 2005gj. 
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The Ha flux continues to decrease with time: on day 174 it is ~3 times fainter than on 
day 37 and on day 204 even 5 times. Finally, the last spectrum (day 389) shows no evidence 
of the typical lines of SNe II; at this epoch, the narrow Ha emission should be attributed 
to the host galaxy. This is consistent with SOS and with the upper limit in the optical 
luminosity that was deduced from the photometry. 

The emission peak of Ha remains at the rest frame wavelength at all phases, exhibiting 
a three-component profile (Figure 7). For the intermediate Ha component we measured a 
FWHM~ 2100 km s"^ at a phase of 42 days, and FWHM~ 3200 km s"^ at 174 days. S07 
pointed out an asymmetry of the line at early times (also evident in our day 42 spectrum), 
likely caused by a blueshifted P-Cygni absorption, which vanishes with time. For this reason 
we can admit that the true unabsorbed profile Ha remained roughly constant during the SN 
evolution. A roughly constant FWHM~ 9100km s~^ is measured for the broad component 
of Ha. 

The physical interpretation of the intermediate and broad components is still a matter 
of debate. S07 and SOS assumed that the intermediate component (v ~ 4000 km s~^) traces 
the kinematics of the SN shock wave, while the broader one is related to the SN ejecta 
(v ~ 6000 km s~^, a value significantly lower than what we obtain, v ~ 9100km s-^). The 
intermediate velocity component was used to c ompute the luminosity expec t ed fro m CSM 



interaction. On the other hand, according to I Chevalier &: Fransson I (119941 . 1200 ll ) and to 



Z ampler i et al.l (120051 ) . the luminosity originating from the reverse shock during ejecta-CSM 
interaction is proportional to the ejecta velocity, i.e. to the width of the broadest Ha compo- 
nent. There is still no consensus on this issue. Because of these ambiguities, one should be 
careful before assigning physical velocities to various regions from just line widths, especially 
for objects with peculiar individual features as are SNelln. 



5. Discussion 



As discussed in §3, the light curve of SN2006gy shows three distinct phases: i) a. very 
broad, exceptionally high luminosity peak (day to ~170), ii) an intermediate phase of slow 
decline (day ~170 to ~237) and Hi) a late phase in which the optical luminosity drops below 
the detection limit and IR emission dominates (day >3S9). As we will show, the first phase 
requires a specific star+CSM configuration. The other two phases have been observed in 
other SNe. 

In the following, we rewind the movie of the event and use the late observations to 
constrain the possible scenario. Starting from the late phases, we discuss the role of dust 
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and ^^Ni in the ejecta, stressing that a very large amount of ^^Ni is not required. We then 
consider the evolution of the SN at intermediate phases and explain that, independently of 
the source that powered the luminosity at peak, interaction dominates between days ~ 170 
and ~ 237. Finally, we discuss the light curve models for the first ~5 months obtained with 
a semi-analytical code (jZampieri et al.ll2003l ). Based on these results, we propose a new 
evolutionary scenario for SN 2006gy. 



5.1. Nickel mass and dust emission 



The late light curve of most SNe is powered by the radioactive decay of ^^Ni into ^^Co 
and ^^Fe via 7 and e"*" deposition. 



Thermonuclear SNe l a eject a large ^^Ni mass (O.IM0 < MssNi < l.lM0. ICappellaro et al 



19971 : iMazzali et al.ll2007l ). but become rapidly transparent to the 7-rays from the radioactive 
decay because of the small ejected mass and the high expansion velocity. As a consequence, 
at t ~ 100 days after explosion the luminosity declines at a rate ~1.5mag (100 d)^^, higher 
than the ^^Co decay input 0.98mag (100 d)~^). A similar behavior was found for most 



type Ib/c SNe fIClocchiatti fc Wheeled 119971 ). 



In the case of H-rich, core-collapse SNe the ejecta remain almost opaque to 7-rays for 
more than a year, and the late-time luminosity decline tracks the radioactive decay. In 
this case, if the date of the explosion is known, the late-time luminosity provides a direct 
estimate of the ejected ^^Ni mass. This spans a wide range of values (0.005 < Mse^i < O.SM©; 
L. Zampieri et al., in preparation), but is typically smaller than in SNe la and Ib/c. 

For SN 2006gy, in the optical bands only upper limits to the luminosity at very late 
phases (411 days) can be obtained. In the near infrared, the K-band detection reported in SOS 
and discussed in the previous section may be suggestive of the presence of low-temperature 
dust emitting in the far IR. This makes a precise estimate of the ejected ^^Ni mass difficult. 
The bolometric luminosity including the emission from dust (plus symbols in Figure 4) imply 
ejected ^^Ni masses up to ~ ISM© for T=800K. However, given the uncertainty on the 
nature of dust and its temperature, a more significant estimate of M(^^Ni) can be obtained 
adopting the bolometric luminosity at earlier epochs, i.e. at ~ 180 days. At this phase the 
relation L = 1.4 • IO'^^Mni • exp(— 1/113.6) ergs~^ provides M^i ~ 3Mq assuming complete 7- 
ray trapping (see also § 5.3). This value is in disagreement with the value obtained by 
Smith fc Mc Crayl (120071 ) with the same relation (MseNi ~ SM©). A possible explanation 



may resi de in a different es t imate of the bolometric luminosity, which is about a factor 3 
higher in ISmith fc Mc Cravl J2OO7I ). 
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Of course, we expect the bolometric flux - and therefore the ^^Ni mass - to increase if 
the IR/longer wavelength emission contribution is taken into account. On the other hand, 
the luminosity decay at phase 170-137 days is much slower than what is expected from ^^Co 
decay. This suggests that an energy source additional to radioactive decay of ^^Ni has to be 
present. Compared to those measured for other SNe, an amount of 3M f^ of ^^Ni may not 
appear unreasonably large (see for example SN 1999as. iDeng et al.ll200l[ ). 



5.2. Evidence of strong, late-time ejecta-CSM interaction 



We mentioned in §3.1 and §4 that at 170-237 days SN2006gy shares several properties 
with SNel997cy, 1999E, 2005gj and 2002ic. Although some authors reg ard some of these 
SNe as thermonuclear expl osions fsee lHa.muv et al.ll2003l for SN2002ic and lPrieto et al.ll2007 
for SN2005gj, but see also lBenetti et al.ll2006l and S. Benetti et al. 2008, in preparation, for 
an alternative scenario) there is unanimous consensus on the fact that interaction dominates 
their emission at late phases. Despite the brighter magnitude at maximum, SN 2006gy has 
luminosity and luminosity decline rate comparable to the SNe mentioned above at 170-237 
days (Figure 4), which is when they also show similar spectra. 

Therefore, it is natural to assume that at this phase ejecta-CSM interaction plays a 
dominant role also in SN 2006gy. Although the low X-ray flux at this phase (cf. §1) might 
appear to be in contradiction with the ejecta/CSM interaction scenario, this may not be a 
problem, because for sufficiently high de nsities {p ~ lO^gcm" ^) the X-rays that are produced 



in the shock are immediately absorbed (ITuratto et al.ll2000l ). 



In the context of interaction, the luminosity L arising from the shock is proportional to 
the progenitor mass loss rate M, to the ejecta velocity Vej and to the unshocked CSM wind 
velocity Vcsm, as follows: L oc Vj^jMV^gj^. Unfortunately, because of the ambiguity in the 
interpretation of emission line profiles (§4), we cannot precisely measure the velocities in the 
different circumstellar regions and thus derive a reliable estimate of the CSM density from 
the observational data. However, the emission lines in SNe 1997cy and 1999E are generally 
broader than in SN2006gy (i.e., their ejecta are probably faster), but their luminosity is 
comparable. In view of the former relation, we expect that the shock wave of SN2006gy 
encounters a higher CSM density at 170-237 days. 
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5.3. A highly energetic supernova impinging on massive gaseous clumps 



The SN evolution during; the first 170 days is explained reasonably well by the scenario 
proposed by ISmith fc Mc Crayl (120071 ). In the shocked-shell diffusion model the supernova 
light is produced by diffusion of thermal energy after the passage of the SN shock wave 
through a shell of lOM© of material, ejected in the decade preceding the explosion. The shell 
is supposed to be initially optically thick and acts as a pseudo- photosphere, so that the long 
duration of the peak and the weakness of X-rays emission are explained in terms of a long 
diffusion time and a very large optical depth. The interaction features typical of type Iln 
SNe are supposed to arise in the observed spectra as soon as the blast wave breaks out of 
the opaque shell into the surrounding, lower-density wind. 



However, in ISmith fc Mc Crayl (120071 ) a number of items have not been considered. A 



f irst inconsi stency concerns the model assumptions. Acco rding to the model of iFalk fc Arnett 



( 1 19731 . 119771 ) which is adopted in lSmith &: Mc Crayl (120071 ) . in order to reproduce the observed 
luminosity rise to maximum, the initial radius of the shocked shell has to be much sna aller 
than the radius at peak luminosity. However, in the model of ISmith fc Mc Cravl (120071) the 
initial and final radius differ only by a factor of 2. Thus the model of iFalk fc ArnettI (119731 ) 
is not applicable: the simple assumption of the existence of a single shell at a large radius 
surrounding the exploding star can not explain the properties of the light curve of SN 2006gy, 
in particular the slow rise to maximum. 



Secondly, in the model of ISmith &: Mc Cray! (120071 ) the important role of ^^Ni is over- 
looked. No attempt has been made to estimate the amount of ^^Ni deposited by the SN and 
to determine its effect on the light curve during the diffusive phas^. 

The third problem concerns recombination, whose effects can not be neglected as soon as 
the decreasing photospheric temperature reaches the gas recombination temperature during 
the post-diffusive phase. 

With these shortcomings in mind, we have developed an alternative, comprehensive 
scenario that attempts to take all these aspects into account. First of all, we divided the 
evolution of SN 2006gy into two distinct phases, before and after maximum luminosity. Each 
phase was modeled independently. The earlier phase (i.e., the rising branch of the light curve) 
was modeled as the explosion of a core-collapse SN originatin g from a compact progen itor. 
For the peak phase we adopted a scenario similar to that of ISmith fc Mc Crayl (120071 ) . in 
which the ejecta impact on very massive (> 6M0, see Table 5) clumps of previously ejected 



The estimate of 8M0 of Ni reported in their paper derives from the extrapolation of the hght curve 
himinosity after day 170 (cfr. §5.1) 
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material and deposit their kinetic energy. Because the density is very high, the energy of the 
shock produced by the ejecta-clump interaction is completely thermalized. A photosphere 
forms, so that the evolution of the shocked clumps can be modeled as if it was another SN 
with very large radius and little ejected ^^Ni. 



A fundamental difference with respect to the model of ISmith fc Mc Cray! (120071 ) is that, 
in our scenario, the true SN explosion is not completely hidden by the circumstellar material 
which is therefore not homogeneously distributed around the star. Rather, it is fragmented 
into big clumps which may be distributed symmetrically with respect to the centre of the 
star. This is motivated by the assumption that the progenitor of SN2006gy may have 
undergone mass-loss episodes similar to those observed in rj Carinae. The rise to maximum 
corresponds to the early emission of the SN ejecta during t he initial phase in which the radius 



is rapidly expanding, similarly to the case of SN1987A flWoosleylll988l ). In our model the 
peak luminosity is sustained by the combined contribution of the early SN explosion and the 
energy from the ejecta-clump interaction. Unfortunately, no early spectra are available to 
verify this claim. The first available spectrum (37 days) already shows signs of interaction, 
mainly in the Ha profile, probably caused by flux arising directly from the interaction, being 
not thermalized by the dense clumps. Therefore we can reasonably assume that at this phase 
the ejecta-clump collision had already started. Another assumption of our model is that the 
impact is instantaneous, i.e. that all material is reached by the ejecta at the same radial 
distance from the star. 



Our semi-analytical code (see lZampieri et al.ll2003l for more details) was used to estimate 



the parameters of the ejected envelope from a simultaneous comparison of the observed and 
computed light curve, photospheric gas velocity and continuum temperature. The radius 
of the star at the explosion, the mass and velocity of the ejecta and the explosion energy 
are fitting parameters, whereas the ejected ^^Ni mass is an input fixed parameter, which 
is based on the late-time light curve. The fitting parameters are estimated by means of 
a minimization procedure for both evolutionary phases (i.e. the SN explosion and the 
ejecta-clumps impact). 

The parameters of the models for each phase are listed in Table 5 and 6. Models of the 
earlier phase {el, e2, e3 and e4) refer to different values of the input parameters M^i and 
Tree, while models of the later phase (ci and c2) refer to different minima. 

Critical parameters for the earlier phase are the initial radius and the mass of ^^Ni. As 
discussed before, the large increase in luminosity in the pre-maximum phase calls for small 
initial radii (< 10^'^ cm), which are not compatible with RSG stars but are consistent with 
BSG or Wolf Rayet stars. The amount of ^^Ni determines the peak luminosity. The adopted 
upper limit is Mm ~ 2M0, considering that ~ 3Mq were estimated at ~ 180 days neglecting 
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the contribution of interaction, which instead is aheady active at that phase, as explained 
above (according to the model the 7-rays trapping is always more than 80% effective at a 
phase 170-237 days, given the large ejecta masses). A minimum ^^Ni yield of 0.75 Mq was 
required to fit the early rise of the light curve, assuming no contribution from interaction 
(i.e., the CSM is supposed to be rarified in the vicinity of the exploding star). 

Table 5 lists the best-fit-parameters for the initial radius of the star, for the ejecta 
mass and for the velocity and SN explosion energy, given the adopted ^^Ni masses and 
recombination temperatures. It should be noted that the radius estimated by the code is 
actually an upper limit: the initial part of the light curve is not very sensitive as long as 
it remains below the reported values. The SN parameters are not exceptionally high for a 
core-collapse SN. For example, the explosion energy is only ~ 3 — 4.5 times larger than that 
of SN 1987A. The explosion energy increases with increasing ^^Ni mass, as one may naively 
expect from the fact that larger amounts of ^^Ni may be synthetized in more energetics 
events. On the other hand, for a constant ^^Ni mass a smaller recombination temperature 
implies an increase in the ejecta velocity and mass, and therefore in the explosion energy. 

The later phase is not powered by ^^Ni alone. The main source of energy is in fact the 
transformation of the kinetic energy of the ejecta into thermal energy and radiation inside 
the dense clumps, which form a photosphere. The duration and shape of the luminosity 
peak depends on the radius and mass of the clumps and on their expansion velocity. The 
parameters listed in Table 6 are the clump radius, mass and velocity, the amount of ^^Ni in the 
clumps, the recombination temperature, the energy released by the ejecta-clumps interaction 
and the diffusion time. The recombination temperature adopted is T = 6500 ± lOOOK, as 
measured from the 37-day spectrum. For both models reported (ci and c2) the energy 
deposited by the ejecta in the CSM is about a factor ~2-30 smaller than the SN explosion 
energy. This value may result naturally, considering that the clumps cover a solid angle not 
larger than 27r as seen from the centre of the star. In the two models the radius of the clump 
is significantly different. For an ejecta velocity of 8000 kms"^, and assuming that the ejecta- 
clump impact occurs at ~30-40 days, the distance of the clumps is ~ lO^^cms"^. Adopting 
a characteristic sound speed of ~ lO'^cm"^, the shock wave produced by the impact takes 
~ 100 days to cross the clump in model cl and ~ 10 days in model c2. On these grounds, 
model c2 seems to be favored, as the optical display of the shocked clumps is fully developed 
by ~ 40 days after explosion. The val ues obtained for th e clurn p distance and mass are 
roughly consistent with those derived by ISmith fc Mc Crayl (120071 ) . 

Our simple model gives a satisfactory fit for both the explosion and collision phase 
(Figure 8). We did not attempt to fit the light curve in the transition phase. The parameters 
that characterize the explosion of SN 2006gy are actually not particularly remarkable. An 
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extra-ordinary amount of ^^Ni in the ejecta is not necessary to fit the hght curve. The 
estimated amount of ^^Ni i s 2 to 6 times larger than that derived f or other well studied, 
bright core-collapse events (ITuratto et al.l 120001 : iMazzali et al.ll2006l ). It should be noted 
that an high amount of ^^Ni is even not to relate to the huge bright ness of the recently 
discove red SN 2008e s , the second most luminous SN known, according to lGezari et al.l (120081 ) 
and to iMiller et al.l (120081 ). SN2006gy was certainly a highly energetic event compared to 
other normal CC-SN e, perhaps comparable to the class of hy pernovae (e.g., SN20 03dh, 
Matheson et"aDl2003l : SN2003id. Ivalenti et allbood : SN 1998bw. [iwamoto et a"l1ll998h . The 
combined mass of the ejecta and of the clumps is ~ 2OM0, indicating an originally very 
massive progenitor, likely much more massive than ~ 30Mq if account is taken of the likely 
large mass loss in the pre-SN stage. Still, these values are significantly smaller than those 
claimed in some of the previously proposed scenarios (> IOOMq). 



There is increasing evidence for the association of bright SNe I In and LBVs (ISalamanca et al 



2002; Smith fc Owocki 2006; Kotak fc Vink 2006; Gal- Yam et al. 2007; Trundle et al. 



20081). 



The properties of these events seem to require that their progenitor stars experienced mass 
loss rates of the order of ~ O.IMqJt'^, which are only compatible with those estimated for 
stars like 77 Carinae. Also in this case strong, LBV-like mass-loss phenomena are required to 
produce massive clumps around the star. Given the radius at explosion derived by the model, 
a star in a LBV or early Wolf-Rayet phase might be good candidates for the progenitor of 
SN2006gy. 



6. Summary and conclusions 

New observational data of SN2006gy allow us to derive constraints on the physical 
processes underlying SN 2006gy. We confirm the luminosity drop in the R band at days 362 
and 394 first reported by S08, and find a similar drop also for the B and I bands at similar 
epochs. The absence of SN features in the spectrum at t ~ 389 days supports this. In all 
bands the light curves exhibit a broad luminosity peak at day ~ 70, followed by a steep 
decline and then by a fiattening at ~170-237 days. At very late phases the SN was detected 
in the K band. This may deal either with dust formation in the ejecta or with IR echoes 
events. However, the uncertainties on the late-phases scenario do not hamper the estimate 
of ^^Ni mass ejected. Based on the bolometric luminosity at day ~180 un upper limit of 
~ 3M0 is derived. At this epoch interaction, rather than radioactive decay, has been proved 
to be the dominant source. 

During the first 3 months the behavior of SN 2006gy can be reproduced as the explosion 
of a compact progenitor star (with explosion energy ~ 4 — 9 ■ lO^^erg, R ~ 6 — 8 ■ lO^^cm). 
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The SN ejecta collide with some previously ejected material (~ 6 — IOMq) distributed in 
highly opaque clumps. The increasing size of the SN ejecta, the relatively large amount of 
^^Ni ejected, the collision with the extended, opaque clumps with a long diffusion time are 
the "ingredients" responsible for the slow increase of the light curve to maximum and for the 
brightness and extension of the peak. The values derived for the mass of t he clumps and their 
radial distance are consistent with those derived for the shocked shell by ISmith fc Mc Cray 
(l2007h . 



The spectra of SN2006gy at ~170-237 days are similar to those of a number of bright, 
interaction-dominated SNe (Figure 6), with which SN2006gy shares remarkable photometric 
similarities (Figure 4). This confirms that at this epoch interaction plays a dominant role 
also in the case of SN2006gy. In the massive-clump scenario, it is likely that interaction 
signatures start to dominate the spectra when the clumps become transparent because they 
recombine. The ejecta encounter regions of progressively lower density with time. At the 
epoch of our last optical observation (~423 days), interaction has probably already ceased. 
This is consistent with the non-detection of Ha in the Keck spectrum of SOS. 

In conclusion, although SN2006gy was very luminous and energetic (~ 3 to 4.5 times 
more energetic than SN1987A), it does not appear to be an extraordinary event. In fact, 
neither the explosion of a supermassive progenitor, nor extremely high Ni-rich ejecta are 
required to explain the observations. 

Unfortunately, the nature of the progenitor star of SN2006gy still remains obscure. 
Nevertheless, to account for such a violent and energetic explosion as well as for the existence 
of an extremely dense CSM around the exploding star it is natural to consider a very massive 
progenitor, likely more massive than 3OM0, which experienced strong mass loss episodes just 
before the explosion. 

LBV-like mass loss rates resulting in highly dense circumstellar shells seem to be com- 
mon near t he bright end of int eracting sup ernovae (e.g., SN20 06jc, iPastorello et al.l 12006 



20081 and S . Bene tti et al. 2008, 



Smith et al.l (120081 ) to explain the 



SN 1997eg, iHoffman et al.l l2007l : SN 2005gj , iTrundle et al. 
in preparation). A LBV-like outburst was also claimed by 
observational data of SN 2006tf. Considering the radius of the progenitor required to explain 
the long luminosity rise to maximum, the progenitor of SN 2006gy was probably a LBV or 
an early Wolf-Rayet star. 
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Table 1. Journal of photometric and spectroscopic observations of SN2006gy. 



UT Date 


Telescope 


Equipment Bands 


Grisms 


Spec, range [A] Resolution 


[A] Pixelscale ["/pix.] 


/no /OK 
Ub/Uy/zo 




DOLORES 




LlX-D, Lix-tx 


3200-9000 


18, 17 


n OK 
U.zo 


06/09/30 


NOT^ 


ALFOSC 


BVRI 


4 


3400-8800 


21 


0.19 


06/10/29 


Ekar 1.82m 


AFOSC 


BVRI 








0.46 


06/12/19 


Ekar 1.82m 


AFOSC 


BVRI 


4 


3500-7500 


24 


0.46 


07/02/10 


NOT 


ALFOSC 


BVRI 


4,5 


3500-9800 


21,'20 


0.19 


07/03/10 


Ekar 1.82m 


AFOSC 


BVRI 








0.46 


07/03/12 


Ekarl.82m 


AFOSC 




2,4 


3400-7600 


38, 24 


0.46 


07/04/13 


Ekar 1.82m 


AFOSC 


VRI 








0.46 


07/09/14 


Ekar 1.82m 


AFOSC 


R 


4 


3400-7600 


24 


0.46 


07/10/05 


TNG^ 


NICS 


JHK' 








0.25 


07/10/17 


TNG^ 


DOLORES 


BRI 








0.25 


08/01/12 


TNG^ 


NICS 


K' 








0.25 



'^Nordic Optical Telescope 
^Telescopic Nazionale Galileo 



Table 2. Main data on the archive images used for the photometric template subtraction. 



UT Date Telescope Equpment Bands Exp. time Seeing ["] Pixelscale ["/pix.] 



91/12/01 JKT'^ 


AGBX 


B 


600 


1 


0.33 


91/12/01 JKT 


AGBX 


R 


300 


1.15 


0.33 


96/01/13 JKT 


AGBX 


I 


360 


1.5 


0.33 


03/02/12 Schmidt Tel. 




V 


300 


3.7 


1.46 



^Jakobus Kapteyn Telescope 
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Table 3. Optical photometry of SN 2006gy. 



UT Date 


JD -2,400,000 Phase [days]'' 


B 


Berr 


V 


Verr 


R 


Rerr 


I 


lerr 


06/09/18 


53996.5 


29.5 










15'' 








06/09/30 


54008.5 


41.5 


16.00 


.08 


15.19 


.22 


14.51 


.10 


14.48 


.07 


06/10/29 


54037.5 


70.5 


15.84 


.06 


14.85 


.13 


14.28 


.04 


14.10 


.06 


06/12/19 


54088.5 


121.5 


16.23 


.06 


15.08 


.15 


14.99 


.04 


14.37 


.06 


07/02/10 


54142.4 


174.5 


17.87 


.08 


16.75 


.22 


16.69 


.10 


15.86 


.07 


07/03/10 


54169.5 


204.5 


18.07 


.06 


17.74 


.15 


16.67 


.04 


16.19 


.06 


07/04/13 


54203.5 


236.5 






17.66 


.15 


16.93 


.04 


16.32 


.06 


07/09/14 


54356.5 


389.5 










>20.30 








07/10/17 


54390.6 


423.5 


>21 








>21.55 




>19.75 





*With respect to JD=2453967.0 
''Prom ATEL 644 (2006) 



Table 4. Near-infrared photometry of SN 2006gy. 



UT Date 


JD -2,400,000 Phase [days]^ J Jerr 


H Herr K 


Kerr 


07/10/05 


54378.5 411 > 17 - 


> 16.5 - 16.00 


.50 


08/01/12 


54477.5 510 


- 16.3 


.50 



^With respect to JD=2453967.0 
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Table 5. Model output parameters of the semi-analytical code for the first evolutionary 

phase. 



model Rstarl 




Mej[M0] 


Vej[km/S] 


MNi[M0] 


Trec[K] Eexpi[105ierg] 


(el) 


8.4 


5.3 


7700 


0.75 


7000 3.8 


{e2) 


5.9 


8.3 


8900 


1.0 


6500 7.9 


(eS) 


8.4 


6.9 


7700 


1.0 


7000 4.9 


{e4) 


8.4 


UA 


7300 


2.0 


7000 9.2 



Table 6. Model output parameters of the semi- analytical code for the second evolutionary 

phase. 



model Re,[-10i2cm] Mei[M0] Vei[km/s] Mm[MQ] T,^,[K] Ei^pflO^ierg] diff. time [days] 



(cl) 1339.9 6.5 1900 0.1 6500 0.3 100 

{c2) 290.6 10.0 3600 0.1 6500 1.6 10 
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Fig. 1. — BVRI absolute light curves of SN2006gy, obtained with the distance and extinction 
reported in t he text. Late phase (>3 00 days) dete c tion l imits are marked with an arrow. 
R data fro r n ISrnith et al.l (120071 ) and ISmith et al.l (120081 ). as well as R and I data from 
Ofek et al.l (120071 ) are also reported. 
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Fig. 2. — Images aquired with NICS at TNG with filter J (left panel) and K' (right panel) 
on October 5th, 2007 (JD 2454378.5). SN 2006gy is still clearly visible near the host galaxy 
nucleus in the K' band image, whereas there is no detectable source at that position in the 
J band frame. 
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Fig. 3. — Comparison between the optical (day 423) and near-infrared (day 411) flux mea- 
sured for SN 2006gy. We also show the expected emission from dust at different temperatures, 
normalized to the K band magnitude. 
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Fig. 4. — Pseudo-bolometr i c ligh t curve of SN2006gy c ompared to t hose of type IIP 
SN1987A Jwhite fc Mali nI Il987h. t ype Iln SN20 05gi dPdeto et all [2007i) . SN 1999E 
jRigon et al.l boosh . SN 1997cv jTuratto et al.l bood : iGermanv et all bood ) and SN 1995G 
(Pastorello et al. 2002), all integrated in the same wavelength range. Red crosses at late 
times include the near-IR contribution due to a possible cold dusty region in SN2006gy 
ejecta, based on the K-band detection and on three possible dust temperatures (T1=800K, 
T2=1000K and T3=1200K, see §3.1). For SN1997cy and SN 1999E the epochs of the as- 
sociated GRB explosions (GRB 970514 and GRB 980919) are adopted as phase reference 
epochs. 
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Fig. 5. — Spectroscopic evolution of SN2006gy from 37 days to 389 days since explosion in 
the host galaxy rest-frame, corrected for extinction assuming E(B— V)=0.56. The spectra at 
phase 174, 204 and 389 were multiphed by a factor 2, 3 and 6 respectively. 
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Fig. 6. — GEL ATO comparison bet ween the spectrum of SN 2006gv at phase ^ 174 days, 
SN1997cy ftop. IXuratto et aLl[200o[ ) and SN 1999E fbottom. [jigon et all [2OO3I ) at similar 
phases. Ahhough the c omparison SNe have broader hnes, (e.g., FWHMHa=12800 kms~^ in 
SN 1997cy according to lTuratto et aLll2000l ). the objects show an overall remarkable similar- 
ity. 
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Fig. 7. — Detail of the Ha profle in the spectrum obtained at NOT on February 10th, 2007. 
The hne is decomposed into three gaussian profiles, having FWHM = 685 (unresolved), 3200 
and 9000 kms'^ 
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Fig. 8. — Best fits of the light curve of SN 2006gy obtained with the semi-analytical model 
(jZampieri et al. I l2003h . showed separately for the rising branch and maximum/post maximum 
phase. The code in the legenda refers to the models summarized in Table 5 and 6. 



